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Abstract 

> 

^-j. ■ We work out finite-dimensional integral formulae for the scalar product of genus 

Q\ ' one states of the group G Chern-Simons theory with insertions of Wilson lines. As- 

suming convergence of the integrals, we show that unitarity of the elliptic Knizhnik- 
Zamolodchikov-Bernard connection with respect to the scalar product of CS states is 
closely related to the Bethe Ansatz for the commuting Hamiltonians building up the 
connection and quantizing the quadratic Hamiltonians of the elliptic Hitchin system. 



1 Introduction 



The present paper continues the program [17][10][19||[20| aimed at analysis of the scalar 
product of states in the Chern-Simons (CS) theory. It extends the considerations of ref. || 
were we treated the SU(2) CS theory on the elliptic curve (times the time-line and with 
insertions of time-like Wilson lines) to the case of a general group G. As in the previous 
papers of the series, the point is to express the formal scalar product of the CS theory, 
given by a functional integral over gauge fields, as a multiple finite-dimensional integral. 
The latter, if convergent for every state, provides the space of CS states W with a Hilbert 
space structure and the holomorphic vector bundle W, obtained by varying the modulus 
r of the elliptic curve and the positions z n of insertions, with a hermitian structure. 

The integral expressions for the scalar product of the CS states are close cousins of the 
contour integral expressions for the conformal blocks of the corresponding WZW conformal 
theory. In the elliptic case, the contour integral representations were recently studied in 
ref. IQ. Our approach elucidates the origin of the complicated expressions which appear 
in such representations for a general group: they are induced by a simple trick, already 
used in [1C], which handles a change of variables in the functional integral. 
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The WZW conformal blocks are holomorphic sections 9 of the bundle W of the CS state 
spaces satisfying in the elliptic case the Knizhnik-Zamolodchikov-Bernard (KZB) equations 

Hi 



(8 T + ±H )9 = 0, (d Zn + ±H n )9 = 0, n = l,...,N. (1) 

Above, H n = H n (r,z) are operators acting on the CS states and ft is a coupling constant. 
The KZB equations may be interpreted as the horizontality equations for a KZB connection 
in bundle W. The consistency of the equations (the flatness of the KZB connection) requires 
that the operators H n , n > 0, commute for fixed (r, z). In fact, H n , n>l, are quantum 
versions of the quadratic classical Hamiltonians h n of the elliptic Hitchin system which 
Poisson-commute. In the case of elliptic curve with no insertions the conformal blocks 
coincide with the (linear combinations of) characters of the integrable representations of 
the affine algebras. Hq in eq. ([!]) is then proportional to the Laplacian on the Cartan 
algebra and the KZB equations reduce to the well known heat equation for the elliptic 
theta functions. In the case of elliptic curve with one insertion, Hq(t) becomes a version of 
the Calogero-Sutherland Hamiltonian Q. For many insertions, operators H n (T,z), n > 1, 
are elliptic versions of the Gaudin Hamiltonians [15]. 

A hermitian structure on a holomorphic vector bundle induces a unique unitary con- 
nection of the 1,0 type. One expects that the scalar product of the CS states induces 
this way the KZB connection. For the elliptic case with no insertions, that was shown 



implicitly in [21| where it was proven that the affine characters form an orthonormal basis 
of the space of CS states. For group SU{2) in the elliptic case with one insertion the 
unitarity of the KZB connection w.r.t. the scalar product of the CS states was proven 
in ||]. In the present paper we study this problem for the general elliptic situation. We 
show that, assuming convergence of the integrals giving the scalar product, the unitarity of 
the elliptic KZB connection follows from a result announced in |l^] which provides a basis 
for the Bethe-Ansatz diagonalization of the commuting operators H n (T,z), n > 0. The 
relation of the integral representations for the genus zero conformal blocks of the WZW 
theory to the Bethe Ansatz was observed in [|l| , see also [Q, |2^, [ll]] . The elliptic SU (2) case 
counterparts of these relations go back to XlXth century works of Hermite on the Lame 
operator, as noticed injj6|. Recently, the integral representations for the elliptic conformal 
blocks were used in |l2| to obtain the Bethe-Ansatz treatment of the elliptic Calogero- 
Sutherland model, an open problem till then. Our work exhibits an intrinsic connection 
between the Bethe-Ansatz and the unitarity of the KZB connection. 

The paper is organized as follows. In Sect. 2, we describe briefly the space of CS states 
in the holomorphic quantization and identify the states on the elliptic curve as vector- 
valued theta- functions. Sect. 3 is devoted to the scalar product of the elliptic CS states. 
As realized in [^] for the case with no insertions, the formal functional integral over the 
gauge fields giving the scalar product may be computed as an iterative Gaussian integral. 



The insertions for general group G are handled by combining the methods of refs. [21, 10 1. 
In Sect. 4, we describe the KZB connection and recall its relations to the Hitchin integrable 
systems. Sect. 5 discusses the unitarity of the KZB connection and finally, in Sect. 6, we 
consider the relations to the Bethe Ansatz. 
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2 Elliptic CS states 



Let us recall the description of the CS states on the elliptic curve T T = C/(Z + tZ) with 
insertion points z n , given in [^, H, with due modifications required by the replacement of 
the group SU (2) by G. For simplicity we assume G to be simple, connected and simply 
connected and will denote by g its Lie algebra. Let A denote the (complex) vector space 
of the 0, 1-components A = Agdz of the g-valued gauge fields. The group Q of complex 
gauge transformations g :T T — > G c acts on A by 

A _> 9A = gAg~ l +gdg-\ (2) 

with d = dzdz- The Chern-Simons states ^ are holomorphic functionals on A. For 
the case with insertions of time-like Wilson lines, they take values in the tensor product 
<g) n V\ n =V_oi the irreducible representation spaces of G of highest weights A n , associated 
to the insertions. States verify the chiral Ward identity 

y(0A) = e fe ^-M) 0n g ( Zn ) (n) y(A) (3) 

where S(g,A) is the action of the Wess-Zumino-Witten model coupled to A, see |§, and 
the subsubscript (n) indicates that the group element acts in the factor V\ n of the tensor 
product space V_. 

Restricting functionals to connections A u = ixudzjTi for u in the Cartan algebra h c 
and T2 = Imr, we can assign to every state a holomorphic map 7 : h c — > V_ related to Vl/ 
by the equation 

V(A U ) = e^NVC^) ® n {^ {Zn - 2n)u/T2 ) {n) l{u) (4) 

where \u\ 2 = (u,u) with the Killing form (-, •) normalized so that \a\ 2 = 2 for long roots 
a (we identify g with its dual). Holomorphic maps 7 corresponding to states satisfy the 
following conditions 

7(u + q y ) = 7(11) for q y in the coroot lattice Q v , (5, a) 

7 ( n + rg v } = e -,iK q \ Tq - + 2u) ^ ( e - 2 — <? v ) (n)7 (n), (5,b) 

= ffin(/*)(n)7(w) for (5,c) 
7(untu) -1 )) = 

®n (^)n7(^) for w G normalizing h, (5,d) 



and 



p 



$>- 2 ™^(e a ) (n) 1 (u + tpV)=0(tn (6) 



for any root a, coweight p v satisfying (p v , a) = 1, u s.t. (u, a) = m + rs with m, s integers, 
p = 1, 2, ... and f — > 0. Here and below e a is the step generator of g c corresponding to root 
a. 
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We shall denote the space of maps 7(14) satisfying the above conditions by W TiZi \. 
Conditions (|5[a) and (||,b) mean that 7 is a vector of theta- functions which due to condition 
(||,c) take values in the zero weight subspace V_ C V and which by condition (||,d) are Weyl 
group covariant. Equations (||) specify their regularity on the hyperplanes (u, a) € Z + tZ 
where a € h are roots of the algebra g. They have been obtained for G = SU(2) in || 
and for general G in |7]] , see also [Q . They follow by demanding the regularity at t = of 
the maps 

t i $ V( h ™A <t ) = V(9™A u+tp v) 

where t 1— > A u / tt = ^(u' + tp v + e a )dz, with u' = u — (m + rs)p v , is a 1-parameter 
holomorphic family of gauge fields, h ms (z) = exp[^((m + fs)z — (m + rs)z)p v ] is a 
multivalued gauge transformation and g ms {z) = h ms {z) exp[— t~ 1 e a ] h~^ s (z) is a univalued 
one. Since («', a) = 0, yl u ' +tp v may be gauge transformed to A u ^ t for f 7^ by a constant 
gauge transformation exp[i _1 eo,]). Gauge fields hms A u i$ lie on codimension one strata in the 
space of gauge fields which cannot be attained by gauge transforming A u J s. Conditions (||) 
assure the global regularity of \& but, due to the properties (||,a,b), they are not independent 
for different m and s. For example, for all simple groups different from SU (2) or SO(2r + l) 
it is enough to take m = s = 0. 



3 Scalar product of CS states 



The scalar product of Chern-Simons states is formally given by the functional integral 

|| $ || 2 = J \*(A)\ 2 e-^I^*^ DAD A* (7) 

We shall perform the above functional integration reducing the scalar product expression 
to a finite-dimensional integral which, if finite, will provide W TtZt \ with a natural structure 
of a Hilbert space. The finiteness has been proven in the special cases (general G with no 



insertions [21] or SU{2) with one insertion ||) and has been conjectured to hold in general 
[18|. The strategy for the calculation of the integral (0) will be as in j2l| and @ so we 
shall be brief discussing in detail only the treatment of the insertions essentially borrowed 
from @. 

3.1 Change of variables 

We shall reparametrize the gauge fields in the functional integral (Q) as 

A = a ~\ (8) 

with g a complex gauge transformation and u in some fundamental domain of the action 
of translations u 1— > u + q y , u 1— > u + rq v on h . (h c /(Q v + rQ v )/Weyl group is the 
space of gauge orbits of semistable gauge fields; ignoring the Weyl group action produces 
an overall factor in the scalar product equal to the order of the Weyl group.) We shall use 
the Iwasawa decomposition of the G c -valued field g : 

g = exp[^ v a e a ] expf^/2] U = n exp[0/2] U = bU (9) 

a>0 
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where <fi takes values in the Cartan algebra h and U in the compact group G. Upon the 
change of variables, field U decouples from the functional integral ([/]) due to the gauge 
invariance leaving us with a WZW-type functional integral over the G c /G-valued fields b 

m 

|| * || 2 = J (*(A U ), ®n(bb*){n)( Z n)V(Au)) 

x e kS(bb*,A u +A:)-£ftr(A* u A u ) ]{ ^ h) 8{m) m d 2 u (1Q) 

where Db = f\ z db(z) is the formal product of G c -invariant measures on G c /G, the delta 
function fixes the remaining freedom in the parametrization @ and the Jacobian of the 
change of variables 

j{u,b) = const. T ~2r e 2hVS(bb*,A u +A:) det ^J Au ) . 

Above, r stands for the rank of G, h v for its dual Coxeter number and Ba u = B + [A u , ■ ]. 
The last determinant was computed in [21]: 

det(B\ u B Au ) = const. r 2 2r e" hV|u - u * |2/r2 |n(n)| 4 
where II is the Kac-Weyl denominator: 

n(u) = e nidT/12 Yl (e™<"' a > - e- m{u > a) )Y[{l - q l ) r JJ(1 - ^ e 27ri<u ' a> ) 

a>0 £>0 a 

with d denoting the dimension of the group and q = e 27riT . The WZW action in the 
parametrization @ takes the form 

S(bb*,A u + A:) = ~J (d<j>M -±J (e^(n- 1 Bn u )%-^, n" 1 ^) 

where 

U U = e -n(z-z)u/T2 ne n(z-z)u/T 2 ^ fa = <f, - n( z - z)(u - U)/T2 . 

Note that 

n u = ex P [^e- w{u ' a){z -~ z)/T2 v a e a ] = exp[]T<e Q ]. 
In terms of the Iwasawa variables, the invariant measure on G c /G is 

db = f[ dtf J] d 2 (e~^/ 2 v a ) 

j=l a>0 

where <j>> = (hj,cj)) are the coordinates of eft w.r.t. an orthonormal basis (hj) of h. Using 
the holomorphic functions j(u) to represent \P and the parametrization (0), we obtain 

( V(Au) , H(bb*) - } (z n ) *{A U ) ) = <M Re(|u|2) ( 7 («) , MKr^n-X, (z n ) 7 («) } . 
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Finally, 

-^I tr ( A l A u) = -irk(u,u)/T 2 . 

With all these ingredients, 

||* || 2 = const. /■( 7 (n),®„(«)- 1 e-^n; 1 ) (n) (^) 7 W) exp[-^g^/(90,50)] 



x exp[-^l /(e^K-^nJ^e^, n" 1 ^) + - ^| 2 ] 



x |n(n)| 4 <5(0(O)) d 2r u J] ^ IT d 2 (e~ <<Mz) ' Q>/2 <(2)) . (11) 

j=l a>0, z 

In order to render the n^-dependent terms in the action of the last functional integral 
quadratic, we shall introduce new variables (rj' a ) a> o defined by 

= dr]' a e a . (12) 

Since n u (z + 1) = n u (z) and n u (z + r) = e~ 2mu n u (z)e 2mu , it follows that 

rfoz + 1) = r/ a {z) ^(z + r) = e- 2 ^r,' a (z) (13) 

and that the change of variables is well defined since d is invertible (for generic u) on 
functions satisfying the periodicity conditions (^). It is easy to see that the change of 
variables (v' Q ) h- > (77^) has a triangular nature: 77^ = + F a {{v'p) p< a ) where < a if 
a — /? is a positive root. Hence, the formal volume element does not change: 

n dv <<M * ),a>/2 <(*)) = n ^ 2 (e~^ (z) ' Q>/2 ^(z)). 

a>0, ^ a>0, 2 



Although the action in the functional integral (|ll[), which was polynomial in variables 
i4 (and their derivatives) becomes quadratic in (derivatives of) r)' a , there persists the n u - 
dependence of the insertions in the functional integral (|TT|). Thus, we have to invert relation 
(12) in order to express n u as a function of rj = J2 ila e a- This will be done generalizing 



a trick of 



la- 

a>0 



Let n = Ce a be the nilpotent subalgebra of g . The corresponding group N C G 

may be mapped into the enveloping algebra U(n) (no completion problem arises if we work 
in highest weight irreducible representations of G c ). U(n) is graded by the positive cone 
Q + in the root lattice Q: 

W(n) = U(n) g 

geQ+ 

according to the eigenvalues of the adjoint action of h, with e ai • • • e am € W(n) ai+ ... +am . 
The map z — > n~ 1 {z) = v(z) may be viewed as taking values in U(n) and it satisfies then 
the twisted periodicity conditions 

u q (z + 1) = u q (z), u q (z + r) = e~ 2m ^u q (z) (14) 
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where v q denotes the Z^(n) <7 -component of v. Relation (|l^ ) may be now rewritten as an 
equation 

dzK 1 = -(dzrf)n- 1 , 

for £Y(n)-valued functions and solved with use of the Green functions of the twisted dg- 
operator. Since {n~ 1 )o = 1, we obtain 



= E & M ■ ■ ■ 9J 1 (d- zV ') • 1 . 



(15) 



K=0 



K times 



The Green function of <9j acting on functions obeying conditions (14) may be easily ex- 
pressed by the Jacobi theta function i?i(z) = J2i(-l) l e 7ri( - l+1 / 2 ^ T+27riz( - l+1 ^ satis f ying 



0i(* + l) = -0i(z), i9i(z + r) = -e 
Explicitly, it is equal to 

^ / 1 (0)i?i(x + z 



(16) 



Px(z) 



7r#i(a:)i?i(z) 



(17) 



with x = (u,q). Hence eq. (15) takes the form 

CO , 

"uH*) = Et-^E / P (u,a 1+ ...+a K ){z -yi)(dzVa 1 )(yi) P {u,a 2 +...+a K )(yi ~ Vi) 
K=0 (ax t ...,atc) 

X {9zVa 2 )(y^) p {u,a K )(yK-i-yK)(dzr)' aK )(y K ) d 2 yi ■ ■ ■ d 2 y K . (18) 

In the 7 (it) matrix element, we may insert the partition of unity \^a n , a n)(Ha n -, «n|) 

a n 

where vectors \^a n ^ a n) corresponding to weights fi an form an orthonormal basis in the 
representation spaces V\ n . This gives 

= Y,Ii e ~ {MZn),f * an) \®n(Va n ,a n \ ® n n^ 1 (z„) 7 ( U )}| 2 



into which we may insert the expressions (18) for n~ l (z n ). As in |1C], the final renormal- 
ization of the functional integral over the Cartan algebra- valued field 4> w iH kill most of the 
terms obtained this way leaving only the ones with \/j, an ,a n ) equal to the highest weight 
vectors |A n ) and the sequences of positive roots (ai, . . . ,a p ) composed uniquely of simple 
roots. We may then write 

( 7 (n),® n («)- 1 e-^n- 1 ) (n) (z n ) 7 W) = ]Je-^^ EE / FKa(r,u,z,y) 



K a 



K n 



x II H(. d ^'a nii )(yn,i) d 2 y n ,i (A I ® n (e anl • • • e a „, K j7(t 

n i=l 



(19) 
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where (A| = ® n (A n |, a = (ai t i, . . . , a 1;Kl , a 2 ,i, ><XN,K N ) = (oti,---,a K ) is a se- 
quence of K = J2 n Kn simple roots satisfying 



A' 



N 



s=l n=l 



(20) 



y = (yx t i, yi,Kx, 2/2,1, , VN,K N ) = (yi, • • • , Vk) is a sequence of K points in T T and 

the kernels i^, iQ , are composed of the Green functions of twisted d 

FK,a{T,U,Z,y) = Y[P(u,a n!l +...+a^ Kn )(Zn - Vua) P(u,a n , 2 + ...+a n , Kn ){ynA - y n< 2) 
n 

X ■ ■ ■ P{u,a n , Kn )(yn,K n -X ~ Vn,K n ) ■ (21) 
contains the terms that will drop under the renormalization of the ^-integral. 
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3.2 Functional integration 



The above use of eq.(15) to express the rj dependence of the insertions reduces the rf 
integral to the form 



A 



II (dz'n'as ) (Vs) (dzVp s ) (vs) exp 



-i(fc+2/i v ) 
2^ 



8=1 



a>0 



a>0, z 



A 



( * ) K 



E II e< ' u(Vs),as) 5 ^ is) 5(Vs - v a{s) ) II detCa+ft,)" 1 (22) 



o-eSif s=i 



a>0 



with <T running over the permutations of K points and with d a = e ^ u ' a )/ 2 de^ u ' a ^ 2 . The 
determinants are well known 



H det^da)- 1 = const. \Ii{u)\- 2 [J |1 - q l 



2r 



a>0 



2=1 

x exp[ l -%- f{d<j),d(t)) - ^-\u-u\ 2 }. 



4tt J > r ' 2r2 

After the reintegration and an easy combinatorial manipulation, see pi], trading the 



sum over root sequences a into sums over permutations (two a's satisfying eq. (20) differ 
necessarily only by a permutation), the scalar product formula (11) becomes 



const, e 



-■KTT2/& 



OO - N A 

J J | (/ ; -'' / TT e -(MZn),*n) TT p 

(=1 ^ n=l s=l 

r 

x exp[-^/(^,50) + ^\u-u\ 2 ] 5(0(0)) 11^' l n WI 



2 A 

£ E F K,aa(T,u,z,ay) (A| <g>„ (e ((ra)nl • • • e ((ra)nifn ) (n) 7(u)) d 2r u ]J d 2 y s (23) 

A o-gS x s=l 

where n = k + h y and a is a fixed sequence of K simple roots satisfying (Bel) . 
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The remaining t/>-integral is of the Gaussian (Coulomb gas) type and may be easily 
performed: 

j e J" eX p[-ff !{d<t>M)] 5(0(0)) U 

= const. r 2 r/2 det'(9t5)-/ 2 e "3^ ^ J 

= const. r -^ 2 e^/6 fl |1 - q l \~ 2r e~^) T,, I P(v)S(v-v')P(y') {u) 
1=1 

where we should take f = — J2n ^n^z n + J2s a s$y a which corresponds in the Coulomb gas 
jargon to external charges at points z n and screening charges at points y s . The Green 
function is 

G(z) = In |tfi(z)| 2 + 2^ " + const - ( 25 ) 

Renormalizing the divergences due to the singularities at coinciding points by point- 
splitting, as explained in [jnj and ||, we finally obtain for the scalar product of CS states 
the following finite-dimensional integral expression: 

|| * \\ 2 = const. r~ r/2 f e^t |w -" |2 |e-^ T '^ (G(r, u, z,y) , 9{u) ) ? d 2r u f[ d 2 y s (26) 

where 6{u) = II(it)7(it), 

+ \ H ZnXn ~ \ H y gas > ( 27 ) 



N K 

. 1 

w = u 

K 

n=l s=l 



S is a multivalued holomorphic functions of t and non-coincident z n 's and y s 's and G a 
holomorphic multivalued map with values in the dual space V_* of V given, respectively, by 

S(r,z,y) = ^ (A n ,A n ')ln^i(z n -z n /) - ^(A n , a 8 ) ln??i(z„ - y s ) 

n<ri n - s 

+ J2( a s,a s ')lnMy s -y s ') (28) 

s<s' 

with x (z) = #i(z)/#i(0) and 

G(t, u,z,y) =J2 J2 F K,aa(T, u, z, ay) (A| ® n (e {(7a)n l ■ ■ ■ e {aa)nJ<n ) {n) . (29) 
K <tGS k 



Using the transformation properties (16), it is straightforward to verify that for H = 
e~* s (G, 9) , 5 n = (0, . . . l . . . , 0) and 5 s defined similarly, 

n 

H(T,u + q v ,z,y) = H(T,u,z,y), 

H{T,u + Tq v ,z,y) = ' TqV +2w) H(t,u, z,y) , 

H(r,u,z + 5 n ,y) = (-1)<^>A H(r, u, z, y) , (30) 

H(T,U,Z + T5 n ,y) = (_l){An,A„)M e -^r<A n ,A n )/ K + 2^( W ,A n > H{T,U,Z,y), 

H(r,u,z,y + 5 S ) = (-l)^^/ K H(r,u,z,y), 
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H(T,U,Z,y + T5 S ) = (-1)^'^)/" e -™T(a s ,a s )/ K -2iri(w, a3 ) H( T ,U,Z,y). 



It is then easy to see that the under-integral expression in (^6() is univalued under the 
u i— > u + (r)q v , z n i— > z n + (t)1, y s i-> y s + (r)l transformations. 

It is useful to compare the above expression for the scalar product of the CS states with 
the genus zero ones obtained in |1C]. The genus zero states are determined by their values 
^(0) = 7 at A = belonging to the G- invariant subspace V_ G of the the tensor product V_ 
of the representation spaces. Adapting the notations to the ones of the present paper, eqs. 
(19.20) of |10| read 

r 2 K 

|| || 2 = const. / I e~ S °^ (G°(z,y) , 7 ) I J[ d 2 y s (31) 

J s=l 

where 

S°(z,y) = ^2{\ n ,\ n >)ln(z n - z n /) - ^(A n , a s ) ln(z n - y s ) 

n<n' n,s 

+ ^2(a s , a s ' ) ln(y s - y s , ) (32) 

s<s' 



and 



with 



G°{z,y) =J2J2 Fk{z,<?v) <A| ® n (e ((ra)nl ■ ■ ■ e (aa)nKn ) (n) (33) 

k tr&S K 



FK_(z,y) tj-K IT 2n-|/n,l Vn,l-y n a 2M,A' n -l~»/n,K„ ' 

n 

The similarity to the genus one case is obvious. 



4 Knizhnik-Zamolodchikov (-Bernard) connection and the 
Hitchin systems 



The trivial bundle with the fiber V_ over the space = C N \ A where A contains vectors 
z with coinciding components carries a 1-parameter family of flat holomorphic connections 
defined by 

V- Zn = d Sn , V Zn =d Zn + l -Hl{z) 



where H®(z) are the Gaudin Hamiltonians 15]: 



d t a t a 

H° n U) = EE T^TT ( 35 ) 



a=\n'+n Zn ' 



with {t a ) forming an orthonormal basis of Lie algebra g. Connection V appeared (im- 



plicitly) for the first time in ref. [23]: it was shown there that the genus zero conformal 
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blocks 7 of the WZW theory satisfy the (Knizhnik-Zamolodchikov) equations V7 = 0. In 
fact, the WZW conformal blocks are horizontal sections of a (generally proper) subbundle 
W° C X% xV G . The fibers W 2 ° of W° may be identified by the assignment — * ^(0) 
with the genus zero CS state spaces. The subbundle W° may be described by giving ex- 
plicit algebraic conditions, depending holomorphically on z_, on the invariant tensors in V_ G 
29[ The KZ connection V preserves the subbundles X% x V G and W° of the trivial 



bundle X% x V . 

The extension of the Knizhnik-Zamolodchikov connection to the genus one case was 
first obtained in ref. ||] and elaborated further in ||, [6|, [14|, p!2[ . We shall use the description 
of the genus one CS states by the theta functions 6(u) = 11(14)7(14). The spaces W T>Z of 
states form a holomorphic bundle over the space of pairs (r, z) with no coincidences in z n 's 
viewed as points in T T . The holomorphic sections of W correspond to holomorphic families 
( T :2.) — > 6( T :Z-> ■ )■ The KZB connection is given by the formulae 

V f = d f , V- Zn = d- Zn , V T = d T + \h q (t,z), V 2 „ = d Zn + ~H n ( T ,z), (36) 
compare to eqs. (||). Explicitly |l^, 12 1, 



N 

H {t,z) = ^A u + ( U2 d x p (u,a)(zn ~ z n >) (e Q ) (n) (e_ Q ) (n/) 



n,n'=l ol 



+ ^ J2M Z n ~ ZW? + P'(z n - Z n ,)) b? w h? , (37) 



3=1 



H n (T,z) = - J2 h ln) 9 u j ~ ( Z) lrP («.«>( Z n ~ Zn '1 ( e ")(n)( e -")(„') 

j=l n'y^n ^ " 



+ T, P (Zn-z n/ ))^ n) hQ (38) 



3=1 

where A u = Ej(^) 2 and p = •&' 1 /'&x. The expressions for H n , n > 1, reduce to the 
ones for the Gaudin Hamiltonians (|35|) in the limit r — > «oo. The operators -ff n , n > 0, 



acting, say, on meromorphic functions of u € h with values in j£ commute, see [14] . Their 
commutation forms part of the conditions assuring the flatness of the KZB connection (^) . 
Although the coefficient functions in H n have poles on the hyperplanes (u, a) E Z + tZ, the 
connection maps holomorphic families of genus one CS states into families with the same 
property, due to the increased regularity @ of the CS states on the singular hyperplanes, 
see 



m 



The commuting operators H® for genus zero or H n at genus one are quantizations 
of classical Poisson-commuting Hamiltonians of the, respectively, genus zero and genus 
one Hitchin integrable systems [22, |3|, 24, ||. Let us briefly recall this relation. Given a 



Riemann surface S and a group G, let A denote the corresponding space of 0, 1-gauge 
fields and Q the group of complex gauge transformations, as in the beginning of Sect. 2. 
Both are (infinite dimensional) complex manifolds and we shall work in the holomorphic 
category. The cotangent bundle T* A is composed of pairs [A, $) where $ is a 1, form on 
S with values in Lie algebra g c . The action (||) of Q on A lifts to the symplectic action 
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(A, $) h-> ( 9 yl, 9 <3? = g<J?g _1 ) on T*.A with the moment map 



M(A, <!>) = dA + AA$ + $ AA 

(the 2- forms on £ with values in g c form the space dual to the Lie algebra of Q). If z is 
a sequence of A~ insertion points in S and A a corresponding sequence of Cartan algebra 
elements, then we may symplectically reduce T*A w.r.t. the S-coadjoint orbit of J2n ^nS Zn 
defining the reduced phase space 



i- A 



where Q z \ is the subgroup of Q fixing J2n ^rA„- The symplectic form on T*A descends 
to (the non-singular part of) V s \ turning it into a finite dimensional (complex) symplectic 
manifold (we shall ignore here the singularities of V z ,\)- For a G c invariant homogeneous 



polynomial p on g c of degree d„, the assignment 



defines a map on T*A with values in the space of dp-differentials (sections of the d p -th sym- 
metric power of the canonical bundle on X). All those vector- valued maps Poisson-commute 
since they depend only on $. They descend to V z ,x giving a system of Poisson-commuting 
maps hp with values in the (finite-dimensional) spaces of meromorphic dp-differentials with 
poles of order < d p at the insertion points z n . Their components form a maximal set of 
classical Hamiltonians in involution turning V Z) \ into an integrable system introduced and 



analyzed in 22] for the case without insertions, see g] for the generalization including 
the insertions. 

Let us specify first the above construction to the genus zero case. In that case, (almost) 
each gauge field A is in the gauge orbit of A = 0, i.e. it is of the form 

A = h^dh 

for h G Q with h determined modulo h i— > h^h with constant /io- Equation M(A, <I>) = 
J2n ^n$z n becomes now 

d( h $) = Y^H^Kh- 1 ^) 5 Zn 

n 

which has a unique solution 

n Z Z n 

provided that the sum of the residues vanishes: 

h(z n )X n h~ 1 (z n ) = 0. 

n 

The group Q z \ is composed of arbitrary gauge transformations g s.t. g(z n ) G , acting 
on h by h i— » hg^ 1 . Above, denotes the subgroups of the G c stabilizing X n under the 
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adjoint action. For the symplectically reduced phase space P® x (the superscript referring 
to the genus zero case) we obtain 

PL = {HNeOA„, 5> n = 0}/G c 

n ' 

where \i n = h(z n )X n h~ 1 (z n ) run through the (co)adjoint orbits 0\ n . The latter are nat- 
urally complex symplectic manifolds and it is not difficult to check that, as a symplectic 
manifold, P® x is the reduction of x 0\ n by the diagonal action of G c . Since for an invariant 
polynomial p, p(&) = p( h &), the corresponding Poisson-commuting Hamiltonians on P® , 



arc 



(dz) d » . 



In particular, for the quadratic polynomial given by the Killing form P2 = ( • , ■), we obtain 
the quadratic meromorphic differential 

= -^Y, Az _ zMz _ Znl) <*) 

= Y. {- -, tr^'Ki + —hi) {dz? 

\ 4tt 2 (z - z n y z- z n J 



where the residues at z — z n are 



h n — 27r 2 E 



%n' Zn 



These are, up to normalization, the classical versions of the Gaudin Hamiltonians of eq. 
(|35|) . The latter may be obtained from /i°'s by replacing the coordinates fj,% = (fi n ,t a ) on 
the coadjoint orbit 0\ n by the generators t" of g acting in the irreducible representation 
V\ n obtained by geometric quantization of 0\ n . 

Similarly at genus one, for £ = T T , (almost) each gauge field is in the gauge orbit of 
the gauge fields A u , i.e. it is of the form 

A = h ^A u = (Kh^BiKh) 

with u £ h c and h u = e nz ( u ~ u )/ T 2 , Hence the gauge fields A may be parametrized by pairs 
(u, h) with the identifications 

(u, h) = (wuw -1 , wh) = (u + q v , h q }h) = (u + rq v , h~ q \h) (39) 

for q v in the coroot lattice Q v and w in the normalizer N(h c ) C G c of h c . Equation 
M(A, <3?) = J2n ^n$z n becomes now 

Qihuh)^ = ^(Kh^XniKhy 1 ^) S Zn . (40) 
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Upon decomposing 

(h u h)(z n )X n (h u h)~ 1 (z n ) = ^2^n a e a + M° - Pn 



with in the Cartan subalgebra h , eq. ( |40|) may be solved, provided that X) n ^n = 0> 
with use of the Green functions of the twisted and untwisted 9-operator on T T : 

n \ a / 

where <3?° is an arbitrary constant in h c . The symplectically reduced phase space P z ,\ 
becomes 

P,,X = { {u,<&°,[±) G T*h c x (xOaJ I £< = 0}/iV(h c ) x (Q v + rQ v ) 

n ' 

with the identifications 
(u,$ ,//) = (umu; -1 , rfuT 1 , w/iur 1 ) = (u + (r)g v ,$°, (/i^ v (z„)^ n /i (rk v(z„))). 

As a symplectic manifold, P^a is a symplectic reduction of T*h c x (xO^) by the group 
7V(h c ) x (Q v + tQ v ). The Hitchin Hamiltonians become 

h p {u,<f> ,fj,)(z) 

= P ( $ " + Wi E ( 11 E P <«,°> ( Z " Z «) ° e ° + ^ " Z n) A ) ) (dz) dp 

which for p = P2 and upon writing h P2 (u, fi)(z) = h P2 (u, /j,)(z) (dz) 2 , reduces to 

n n,ri ^ o 

x P <«,a)(2 - Zn)P-{u,a}(z ~ z n>) + p(z ~ Z n ) p(z - Z n >) H° n ,) 

1 

n 

with the residues at z n 

and the holomorphic (actually z independent) piece 

h Q (u,<t>°, E ) = (^°,«& )-^E ( 2 ^E^,«>(^-^'K^ Q 



1 E^'O 21 ~ Z n){Pn, /%) + E ^ ~ ^n) ^n(^, ^) + h (u, (f) , fl) . 



i \ a 
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+ (p(z n - z n ,f + p'(z n - z n >)) (fi° n , . 



Hamiltonians h n , n > are the classical versions of the elliptic Gaudin Hamiltonians H n 
of eqs. @) and (|3f), see |2||5|. 

The SU(n) elliptic Hitchin system corresponding to one insertion has unexpectedly 
appeared recently in the description of the low energy sector of supersymmetric Yang- 
Mills theories 111. 



5 Unitarity of the KZB connection 



One of the essential features of the structures discussed above should be the compatibility 
of the KZ and KZB connections with the scalar product of CS states. The integrals in eq. 
( |3~l| ) have been conjectured in ]l7j (and proved in many cases) to converge precisely for 
invariant tensors 7 € W% and to equip bundle W° with the hermitian structure preserved 
by the KZ connection. The latter condition means that for all local holomorphic sections 
z i — > 7(V) of bundles W , corresponding to holomorphic families z 1— > ^(z) of CS states, 

d Zn \\*f= (*, V^*). (41) 

In order to see why one should expect such a relation, it will be convenient to express the 
scalar product integral in the language of differential forms, following refs. 
oo(y) = - y~ 1 dy. Introduce the V/^-valued K forms 

^°K,g{^y) = u{zi - 2/1,1) Au^yij - j/i j2 ) A . . . Au(y ljKl -l ~ Vl,^) 

A 

/\u(z N - j/at.i) A u(y N> i - y N>2 ) A ... A u(yN,K N -i ~ Vn,k n ) 

x (A| ® n (e anl ... e anK j (n) 

and 



27, 



Let 



n°(z,y) = J2 E (-l) H n K ,aa(z,ay). 

K aeS K 



We may rewrite the scalar product formula (31) as 



const. 



-S° 



Y z 



where Y z stands for the space of y's with y s 7 s not coinciding among themselves and with 
z n 's. We use the conventions that \n°\ 2 = (-l) K ( I< - 1 )/ 2 (^) K n° Af2° and that the integral 
of the forms of degree lower than the dimension of the cycle is zero. Assuming a sufficiently 
strong convergence of the integrals, we may enter with the holomorphic exterior derivative 
under the integral so that 

2 



d II ^ || 2 = const. / d 



L -s° 
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. f9(e-^>» )7 }) A e^>0 i7) 



const 

where the 3-operator on the l.h.s. acts on the z-variables and the one under the integral 
may be taken as acting on both z- and y- variables. The U_*-valued holomorphic multivalued 
forms Q° have two basic properties: 

N 

dfl° = and (dS°) A Q° + ^ dz n A H° n 0° = (42) 

71=1 

with the contragradient action of the Gaudin Hamiltonians on the V/*-valued form Q°. The 
first relation is trivial. The second, more involved one, has been proven in [p8fl . Using eqs. 
(|42|), we obtain 

d ||* || 2 = const, f |e-^°| 2 5>z„A(ft°, (d Zn + ±H n )-f) A e"^ <n°, 7 > 

which implies the relation (f4l|). 

The above analysis has its counterpart for the elliptic case. We again conjecture that 
the integrals in eq. (p6|) converge for 9 corresponding to CS states (in fact exactly when 
the regularity conditions (||) are fulfilled) and that the resulting hermitian structure on 
the bundle W renders the KZB connection unitary, see || for the proof of this conjecture 
for the case of G = SU (2) and one insertion. In order to substantiate the conjecture for 
the case of general G and arbitrary insertions, let us rewrite the genus one scalar product 
integral (^) in the language of differential forms. Define 

LOq(z) = P {Utq) (z)dz - ± d x P {U)q) {z) dr , 



see eq. (17). It is easy to see that div q = 0, if d differentiates only the variables (r, z). Set 

^K,a{T,u,z,y) = u)p lx {z\ — yi : i) A ujf3 12 (yi : i — yxp) A • • • A wp liKi (yi iKl -i - Vx,Kx) 

A ' 

Aw Av,i i z N ~ VN,l) A Vp Nj2 {yN,l ~ VN,2) A ... A U0 NjKn {yN,K N -l ~ VN,K N ) 

X (A I <g> n (e anil ... e anKn ) (ri) (43) 

with p n j = and 

i'=i 

ft(T,u,z,y) = (-^ a ^K,aa{T,u,z,ay) . 

K aeS K 

We still have to dress with the du 3 differentials. The convenient way to do this is to 
define the form 

ft = e~^ s (dw 1 A ... A dw r A (0,0) + ^~Y,dw l A . . .dr . . . A dw r d u , {n, 9)) (44) 
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where w J = (h 3 ,w) with w given by eq. (ETj) and, as before, 6 = II7. It is easy to see that 



eq. (26) may be rewritten as 



t n2 -r/2 / ^\w-w\ 2 

^ || = const. r 2 ' / e^2 1 



2 j2r 
d U 



= const. r~ r/2 / e^ 1 '^ 1 lOp . (45) 



UY T , £ 

This holds since the terms with dr differentials do not contribute to the integral over the 
cycle UY TjZ composed of points (u, y) with identifications u = u + (r)q v , y s = y s + (t)1 
which form the (r+.ff )-dimensional torus T^ +K (with coincidences of y s 's among themselves 
and with z n 's removed). The gain from adding the terms with c^t's is that the form ft 
transforms under the maps !n->«4- (t)<7 v , y s 1— > y s + (r)l (and z n \— > z n + (r)l) exactly 
as the functions H, see eq. (|30|), and, as the result, the form under the last integral in eq. 
( |45|) is well defined on UY T:Z . 

We would like to compute the holomorphic differential of the norms squared || ^ || 2 of 
a holomorphic family of states. Due to the fact that the integration cycles UY TZ depend 
nontrivially on the differentiation variables, this requires a more subtle consideration than 
in the genus zero case where we could ignore this problem (assuming enough convergence 
at the coinciding points). Ignoring again the coinciding points, the geometric setup is as 
follows. We are given a bundle w : M — > Ai of multidimensional tori over the space A4 of 
pairs (r, z). Let us forget for a moment the complex structures on these spaces considering 
them as real manifolds. On we are given a differential form r\ of degree equal to the 
dimension of the fiber (we shall not need an obvious generalization of the story to forms 
of arbitrary degree) and we are studying a function / on the base obtained by fiber-wise 
integration of rj 



f(n) = J i] . 



In local trivialization of the bundle, only the components of rj with differentials along the 
fiber contribute to /. If there are no convergence problems, then obviously 



df(n) = J d ± r j = J d V (46) 
ro -i(W) w-HW) 

where the differential d 1 - in the transverse directions is defined using a trivialization and 
the second equality follows since Sri = (d-d 1 -)] is a closed form when restricted to the 
fibers so that its fiber- wise integral vanishes (again assuming no convergence problems). 

Let us first see how this argument may be used in the simplest case with no insertions 
where we have 

II ,T, Il2 4. -r/2 f %0r\u-u\ 2 , Q , n,2 ,2r . -r/2 f %£r\u-u\ 2 , OI 2 

|| W || = const. t 2 ' / e 2 \"\ u )\ d u = const. r 2 / e 2 |S2| 
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where 9(u + q v ) = 9(u) and 9(u + rq v ) = e mK (i\ rq v +2u) @( u ^ j. e , Q [ s an r-dimensional 
theta-function of degree 2k, and 

r 

ft = e du 1 A . . . A du r + ^— V d u3 du 1 A . . . dr . . . A dif 

is a holomorphic r-form which under u *— > u + (r)g v transforms the same way as the 
theta-function 9. Note the relations: 

dv? aq = ^dTAd uJ n, J2 duj Ad ^ n = 2L dTAA ^ n - ( 47 ) 



The (2r)-form r\ = |fi| is well defined on the bundle M of tori obtained by 

varying with r. Applying the above geometrical considerations to the case at hand, we 
obtain 

d || * || 2 = const. r 2 - r/2 / (d+ £(dr - df)A) (e^ 1 ^' 2 |ft| 2 ) . 

Clearly, only the holomorphic exterior derivative 5 will contribute under the integral to 
d || \I/ || 2 . Explicit differentiation gives: 

d || ^ || 2 = const. r" r/2 | e ^t^' 2 - €Z | 2 dr A + ^ - v?) dv? A 

T 

+ {i-drA +drAd T + ^tfu* A0„* jo) A ?2 

i 

= const. r~ r/2 J e ^\ u - u \ 2 {^\ u -u\ 2 drA + ±- dr A £(u>' - v?) d uJ 

n~ r J 

T 

+ ^- c?t A + dr A d T + ^ dr A A U J ft J A H (48) 



where to obtain the last equality we have used the relations (47). The latter expression 
may be simplified if we notice that 

^drA £, 8 u3 ( e ^ lu - al2 ( {8 uJ + if(v? - &)) fiAfi)= e ^ u ~^ 



f^\u - u\ 2 dr A +^dr Aj2j{u J -u J )d uJ + ^ dr A A u + g^dr A JfJJ A fi (49) 



with e ^ ^(5 u j + ff^M- 7 — u- 7 ) J fi A fi being a (2r)-form well defined on M. It follows 
that the fiber-wise integral of both sides of eq. (|49|) vanishes and, consequently, that 



d T || * || 2 = const. T~ r/2 j e^ |u -" |2 ((d r + J-A u )n)An 

n~ r 

T 

= const. r- r/2 J e^t'"""' 2 9{u) (d T + £A K j 9{u) d 2r u = V T ^) 
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which proves the unitarity of the KZB connection for the case with no insertions in a 
way that maybe was not the simplest (two straightforward integrations by parts would 
do) but which has the virtue of generalizing to the case with arbitrary insertions (modulo 
convergence problems). 



To treat the general case, we shall apply the formula (|46|) to the differential form rj 



I in — in\ 2 



e 



\w—wr —-S i 



with f2 given by eq. (^4|) and satisfying the generalization of relations 

( f47| ) with du J replaced by dw 3 . As before, we obtain, denoting by & the holomorphic 
exterior derivative in all but u directions, 

d||*|| 2 =const.T 2 - r ' /2 J (d+-^dTA)(e^ lw -™ l2 \Sl\ 2 ^ 

= const. r 2 ' / e^ 1 1 II j-2\V) - w\ dr A + — 2_^{ur - w^dw 3 A + ^drA 

UY T , £ V ~ 3 

+& + e duj A On) a n 

= const. T 2 ~ r/2 J e^'^^lw-wfdrA+^drA^-^d^ + ^ dr A 



A large part of the right hand side vanishes by integration by parts with use of the relation 
generalizing (ff9j): 

| |io - w| 2 dr A +^dr A J2j( wj ~ w j )d uJ + ^dr A A u + fi^dr A \ Afi (50) 
and we obtain 

d || 9 || 2 = const. r 2 " r/2 / e^ P ((5' + ^rfrA A u - ^Ai^ n A9 Uj 

= const. r~ r/2 | e ^ w ~^ f(d' + ±^ dr A A u - ± E A£ cte n A e> u , 

+ !E«i4A5 % .)e4 s (fl,9)] A e-^ s (fi,e) d 2r u. (51) 

The d u j -terms may be transformed with use of relation 

K n 

E E #M " Vn,i) A = ^dr A d u ,n, (52) 



n i=l 
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in the notations of eq. (|43| ) and with y n g = z n , which follows easily form the definitions 
of the forms Q and ui q . Eq. (|52|) may be rewritten as 



(53) 



J2Kdz n -^a° s dy s \ AO. = ^drAd^Q - ^dz n A/»J J2, 



with the contragradient action of /i J on the V*-valued form $7. The last relation, upon 

(n) 

substitution to eq. (|5l|), yields 



3 II ^ II = const, r. 



-r/2 



e £|>-™| 2 



e k 



UY r>: 



(d'(n,6) - ±dS A (Q, 9) 



+ J-dTAA u (n,e) - ^dTAY,dvj((d^n,e)) 

j 

+ l - J2 dz n a d uJ ( ( h{ n n , e ) )) a n d 2 ^ 



const. tv 



-r/2 



e 2T 2 



— - ds a si 



UY r , z 



The crucial result are the following equalities: 

#0 = and (95) A ft + tir A tf ^ + X! dZn A = ( 54 ) 

n 

with the contragradient action of H^s. The first of these equalities is a straightforward 
consequence of the closedness of the forms co q from which Q is built. The second, more 
technical one, has been announced in |l2j (as Prop. 9). Using these relations, we finally 
obtain 



d || ^ || 2 = const, tv 



-r/2 



$%\w-w\ 2 



e « 5 



l) K (0 , d'0 + -drAH 



UY T „ 



+ ~ ^<iz n AH n 6)) A ft d 2r u = V*) 



which proves the unitarity of the KZ connection w.r.t. the scalar product pq ) modulo the 
control of convergence of the integrals. 



6 Bethe Ansatz 



The basic algebraic relations (|42| ) responsible for the unitarity of the KZ connection turn out 
to be a disguised (and compact) form of the Bethe Ansatz solution of the eigenvalue problem 
for the Gaudin Hamiltonians (|35[) . This was remarked in [Q] in the context of contour 
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integral representations for the solutions of the KZ equations and developed further in 
26] and in |11] where a relation between the Bethe Ansatz and the Wakimoto realization 
of the highest weight modules of Kac-Moody algebras was explained. 

The elementary fact is that if, for fixed z, a configuration y of K non-coincident points 
in the plane satisfies the equations 

dy s S°(z,y) = 0, s = l,...,K, (55) 

then the second relation of ([42] ) reduces to 

J2dz n A (d Zn S°(z,y) + H° n (z)) n°(z,y) = 0. 

n 

The last equation gives the Bethe Ansatz solution for the common eigenvectors of the 
operators H®(z) acting by the contragradient representation in V_*: 

(d Zn S°(z,y) + H° n (z)) G°(z,y) = (56) 



in the notations of eqs. (32.34). The conditions (55) required for eq. (|5q ) to hold have the 
explicit form 

1 1 

^2(X n ,a s ) — - — = ^2(a s ,, a s ) — . (57) 

n z n Us s ,, g Us' Us 

The above genus zero story has its elliptic counterpart, as remarked in || for the case 
of G = SU(2), see also ||, and in [^] for general simple groups. For fixed r, z and for y 
satisfying the equations 

d ys S(T,z,y) = 0, s = l,...,K, (58) 

or explicitly, with p = n?j 

^2{X n ,a s ) p(z n - y s ) = ^2(a s i,a s )p(y s > -y s ), (59) 

n s'^s 

the second relation of ( |54] ) reduces to 

dr A (d T S(T,z.y) + H (r,zfj + ^dz n A (a 2n 5(r,z,y) + H n (r,zj) ) 

x n(T,u,z,y) = 0. (60) 

Operators H u (t,z), n > 0, act in the space of meromorphic functions / of variable u E h c 
taking values in V^, obeying the periodicity conditions 

f(u + q V ) = /(«), (61) 



with poles possible on the hyperplanes (u, a) € Z + rZ . Eq. (|60D gives the Bethe Ansatz 
solutions for the common eigenvectors of HnS in that space 



(d T S(T,z,y) + Hq(t,z)) G(T,u,z,y) = 0, 
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(d Zn S{T,z,y) + H n (T,zf) G(r,u,z,y) = 
in the notations of eqs. ( ^I|j29| ). 

A slightly modified version of the above argument allows to extend this construction 
and also to find Bethe Ansatz eigenvectors of H^s acting in the space of meromorphic 
Y_Q-valued functions satisfying the twisted boundary conditions 

f(u + q V ) =e«'« V >/(u) (62) 

with £ £ h c . Denoting 

Q = e S + ^lel'r - T,& X n) z n + E^« a ->f« 

ra s 

and combining the second equality of d53) with the relation (53), we obtain 



(dS) A£l + dr A HqQ, + ^dz n A H n VL = (63) 

n 

which results in the eigenvalue equations for the twisted-periodic Bethe Ansatz eigenvec- 
tors: 

(d T S(T,z,y) + H (t,z)) e^G(T,u,z,y) = 0, 
d z J(r,z,y) + H n (r,z)) e^G(r,u,z,y) = 



holding provided that d Vs S = 0, i.e. that 

^(A n , a s ) p(z n - y s ) + (£, a s ) = ^ (a s >, a s ) p{y s , - y s ) . (64) 

n s'^s 

In particular, for £ in the weight lattice P, we obtain more periodic Bethe Ansatz eigen- 
vectors. 

As explained in refs. Q or [12], in the special case of G = SU(n) and one insertion of 
the nm-fold symmetric power of the fundamental representation, operator Hq acting on the 
V_Q-valued functions reduces to the elliptic Calogero-Sutherland-Moser multi-body operator 
and the above techniques were used in 13 1 to obtain its Bethe Ansatz diagonalization. 



The above case is a quantization of the very same elliptic Hitchin system which appeared 
in the effective low energy description of supersymmetric Yang- Mills theories Q. For 
the SU(2) case, one recovers this way Q the classical Hermite's results about the Lame 
operator [EH]. 
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